Fast and robust spin manipulation in a quantum dot by electric fields 
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We apply an invariant-based inverse engineering method to control by time-dependent electric 
fields electron spin dynamics in a quantum dot with spin-orbit coupling in a weak magnetic field. 
The designed electric fields provide a shortcut to adiabatic processes that flips the spin rapidly, thus 
avoiding decoherence effects. This approach, being robust with respect to the device-dependent 
noise, can open new possibilities for the spin-based quantum information processing. 

PACS numbers: 73.63. Kv, 72.25. Dc, 72.25.Pn 



Introduction.- Coherent spin manipulation in quantum 
dots (QDs) is the key element in the state-of-the-art 



spintronics and solid-state quantum information jiol 11 1 . 
Accurate spin manipulation can be achieved by several 
techniques. One of them is the conventional electron 
spin resonance (ESR) induced by magnetic field oscil- 
lating at the Zeeman transition frequency [l[. A more 
robust technique is the spin manipulation with geometric 
Berry phases during the adiabatic motion 0, E| . Nowa- 
days, there is also a growing interest in the electric con- 
trol of spin using spin-orbit (SO) coupling [ljj]. It has 
been applied to high-fidelity spin manipulation on the 
100 ns time scale [7H9(. This highly efficient all-electrical 
method has other advantages. For example, it is easy to 
generate time-dependent electric fields on the nanoscale 
by adding local electrodes and produce spin manipulation 
by making them Zeeman-resonant 0]. As a result, Rabi 
spin oscillations appear at a frequency much smaller than 
the Zeeman frequency making the flip relatively slow and 
prone to decoherence. We shall propose here another 
all-electrical technique to flip spin with high fidelity via 
"shortcuts to adiabaticity" , in a time that can be much 
shorter than any decoherence time. 

Recently, several shortcuts to adiabaticity have been 
put forward to speed up the adiabatic passage of quan- 
tum systems, and achieve a robust and fast adiabatic-likc 
control [13l-|26|. The transitionless or counter-diabatic 
control al gori thms proposed by Demirplak, Rice fl~3| . 
and Berry [14| , are designed to add supplementary time- 
dependent interactions that cancel the diabatic couplings 
of a reference process. The system then follows exactly 
the adiabatic trajectory of the original unperturbed pro- 
cess, in principle in an arbitrarily short time. Tran- 
sitionless quantum drivings have been implemented in 
two- level systems: spins in a magnetic field 14|, atoms 
15 1, and Bosc-Einstcin condensates in optical lattices 
16( . A different shortcut is provided by inverse engi- 
neering the transient Hamiltonian 171 uM using Lewis- 
Riesenfeld invariants 1271. This method has been used for 
time-dependent tra ps Il7|-l 2ll| , atomic transport 22( , and 
other applications [23|, |24(. Although these two meth- 
ods are potentially equivalent [HI , their implementations 
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FIG. 1: (Color online) Schematic diagram of spin dynamics 
of electron in a QD in the presence of electric fields £i(t), and 
perpendicular magnetic field Bo || z— axis. 



and results can be quite different. Here we choose the 
invariant-based inverse engineering approach, since it is 
better suited than the transitionless driving to be pro- 
duced by the desired all-electrical means. 

Model.- We consider the electric control of electron spin 
in a QD formed in the x-y plane of a two-dimensional 
electron gas confined in the ^-direction by the coordinate- 
dependent material composition, under a weak magnetic 
field Bo || z, as shown in Fig. [T] Here the total Hamil- 
tonian H of the electron interacting with the external 
electric field 8 (t) = —dA/cdt is H = Hq + H so + H hlt , 
with [H 



#o = — - + U(x,y) + — o z 

2m 2 

H so = (-ao-y + f3a z )p x + aa x p y , 



PI 



int 



--A(t).v, 



(1) 
(2) 
(3) 



where m is the electron effective mass and o~i (i = x, y, z) 
are the Pauli matrices. Hq represents the kinetic energy, 
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the potential U(x,y), and the Zccman A 2 = g/ig-Boj 
where /is is the Bohr magneton, and g is the Landc 
factor. The eigenfunctions of H are ipj(x,y) \a), where 
|cr) = | ± 1) is the eigenstate of cr 2 , and the spectrum is 
given by Ej ± A 2 /2, where Ej are the orbital eigenener- 
gies in the confinement potential. 

The SO coupling is the sum of structure-related 
Rashba (a) and bulk-originated Dresselhaus (/3) terms 
for [110] growth axis. The vector potential A(i) is in the 
(x, j/)— plane and corresponding spin-dependent velocity 
operators are 



«x = t l H o + H so ,x] = p x /m + pa z - c 
v y = - [H + H SOl y] = p y /m + aa x . 



We focus on the doublet *i = Vi |1),*2 = J -l ) , 
include the higher orbitals by Lowdin partition |28l . |29| , 
and reduce the full Hamiltonian into an effective 2x2 
one [301, 



H 



cff 



9VB 



Z X + iY 
X -iY -Z 



(6) 



where X = B 2 (l + Y = (a//3)(l + Z = 

Bq + (1 + Cx)Bi, with the effects of higher states char- 
acterized by £e and £ y , and the components of the de- 
signed electric field are renormalized by the factors of 
1/(1 + The effective magnetic fields are expressed 
with the SO coupling parameters as B\ = —2e(3A x /cg[iB, 
and B 2 = —2eaA y /cg[iB- The resulting electric fields 
are: 



£x(t) 



gjxs dBi 
2ef3 dt 



£ y (t) 



gUB dB 2 
2ea dt 



(7) 



In practice, some slowly varying electric fields can be ap- 
plied to drive the state from to ^2 adiabatically along 
an instantaneous eigenstate of the Hamiltonian in Eq. 
© ■ To accelerate the driving using the transitionless al- 
gorithm, counter-diabatic fields should be provided [25| . 
However, the common dependence of Y and Z on B\ pre- 
cludes the implementation of the fast driving terms only 
by electric fields. In contrast, invariant-based inverse en- 
gineering naturally leads to an all-electrical driving. 

Dynamical invariant and spin-flip example.- We shall 
design the time dependence of the external electric fields 
to guarantee the state transfer in some fixed time tf 
by using the dynamical 2x2 invariant I(t) satisfying 
the condition dl(t)/dt = ihdl(t)/dt - [H cff (t), I(t)] = 0. 
Parametrizing the Bloch sphere (Fig. [T]), by the angles 9 
and ip, we construct yet unknown orthogonal eigenstates 
|X±(*)> of J(t) as 



\x + (t)) 



cos — e 
sin — 



\x-(t)) 



sm — 

o 2 

cos -e 
2 



-lip 



(8) 



They satisfy I(t)\x±(t)) = X±\x±(t)). Introducing A± 
±gfiBB c /2, we construct the invariant as [25 1 



I{t) 



B r 



cost 



cos< 



(9) 



where B c is an arbitrary constant magnetic field to keep 
I(t) with units of energy. According to the Lewis- 
Riesenfeld theory, the solution of the Schrodinger equa- 
tion, ihdt^ = H cif (t)^, is a superposition of orthonormal 
"dynamical modes", tf(t) = E„ C n e ta "^ \ Xn (t)} ^3, 
where C n are time- independent amplitudes and a n (t) is 
Lewis-Riesenfeld phase, 



a n (t) 



1 



* d 

(Xn(t')\ih— 



H ca (t')\ Xn (t'))dt'. (10) 



From the invariant condition, dl{t)/dt = 0, the angles 6 
and if are related to X,Y and Z by auxiliary equations 



r](X sin ip — Y cos ip), (11) 
■q(X cosipcote + Y smipcotO - Z), (12) 



where 77 = gpLs/h- Since X is a function of B 2 , while Y 
and Z are functions of Bi, once 9 and ip are fixed, Eqs. 
(fTTj) and (fT2|) give the effective magnetic fields 



B, 



Bo 



—(38 cot 8 cosip + (3(ip + r]B ) sin!/? 

77(1 + £r)( a c °l @ — ft s'mip) 
ad cot 9 sin ip + a(ip + i]B ) cos tp — 
77(1 + £, y )(a cot 8-/3 sin ip) 



(13) 
(14) 



from which the electric fields are calculated using Eq. 
J?]). During the spin- flip process, there exist some time 
instants t — t s which satisfy 



a cot 8(t s ) = j3 sin ip(t s ), 



(15) 



and make the denominators of B\ and B 2 zero. To get 
rid of such singularities we impose the conditions 

psmip(t s )[tp(t s ) + riBo-{p/a)8(t s )cosip(t s )^ =0, (16) 
acosip(t s ) ip{t s ) + r)B -(f3/a)6{t s ) cos ip(t a ) =0, (17) 

which make the numerators of B\ and B 2 zero simulta- 
neously. In the following example, we will show how this 
works. 

In general, the eigenstates of the invariant are not the 
same as the instantaneous eigenstates of the Hamilto- 
nian, since I(t) and H (t) do not commute. If we im- 
pose for 6 at t = and tf that 

6>(0) = 0, 6{tf) = tf, 0(0) = 0, 6{tf) = 0, (18) 

then [# cff (0),/(0)] = and [H cS (t f ), I(t f )] = 0, which 
guarantees common eigenstates at initial and final times. 
Moreover the state obeying Eq. (TT5|) will flip from 
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FIG. 2: (Color online) Dependence of the maximum of applied 
magnetic field i?Q lax (solid blue) on the time tf for a third 
order polynomial ansatz for 6 and <p, with the parameters: 
ha = 2 x 10" 6 meV-cm, /3 = a/2, g = -0.44 for GaAs. 
Bo = 0.005 T (dashed red) corresponds to A z — 1.5 mK. 

I^i) at t = to 1*2) at t = tf, up to phase factors, 
along the eigenstate \x+(t)). To design the trajectory 
at intermediate times we assume the polynomial ansatz 
9 = Ylj=o a j^ , where the aj can be fixed by solving the 
system implied by Eq. ([T8]). This leads to 0(t f /2) = tt/2, 
so cot 9 covers the whole (—00, 00) range passing through 
zero at t = tf/2. This may lead to one or several times 
satisfying Eq. (fT5j) . as we will see below in more detail. 

To determine fully, we also need the trajectory 

for tp. As the initial and final states are the poles of the 
Bloch sphere, the phase <p is not well defined there. We 
may nevertheless specify how the trajectory approaches 
them, and impose limits from the right at t = 0, and 
from the left at t = tf, for example, 

^(0+)=7r/2, tp(tj)= w/2. (19) 

These conditions are not sufficient though, since we still 
have to deal with the singularities and their cancellation. 
As cot#(t//2) = 0, we may satisfy Eq. (|15|) and im- 
pose zeros of the denominators for B1.2 at t s = tf/2, if 
smtp(tf /2) = 0. Imposing the two conditions 

<p(t f /2) = 0, (20) 
<p(t f /2) = (P/a)d(t f /2)-TiB , (21) 

to satisfy Eqs. (fll)f and (fTTj) at t 3 = tf/2, we cancel the 
singularity there. With the conditions in Eqs. (|r9|) - (|2T1) . 
we solve the third order polynomial ansatz <p = Ylj=o ^jV 
to determine <p(t). 

Explicit calculations demonstrate that for the third- 
order polynomial ansatz and boundary conditions im- 
posed here, there is only one (removable) singularity at 
t s = tf/2 when the field Bq is smaller than certain up- 
per limit £?™ ax , shown in Fig. [2] as a function of tf. 
For B > B™ ax , more solutions of (|T5|) appear [B and 
siny(i) are coupled by Eq. (f2"Tj)]. which cannot be can- 
celed with the third order polynomial. To satisfy Eqs. 
P6p and (|17|) at more than one zero of the denomina- 
tors of B12, one may set higher order polynomials for ip 
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FIG. 3: (Color online) For tf — 1 ns: (a) Polynomial ansatzes 
of auxiliary angles 8 = X^=o ( s °lid blue), ip — 5Zj=o 
with Bo = 0.01 T (dashed red) and B = 0.07 T (dotted 
black). The designed electric fields £ x (solid blue) and £ v 
(dashed red) by which spin flip can be realized for Bo = 0.01 
T (b) and Bo = 0.07 T (c). Other parameters are the same 
as those in Fig. [2] For simplicity, we put here ^ x = £ y = to 
skip the trivial dependence on these factors. 



and further conditions. This increases the bound i?Q nax , 
but also complicates the driving fields. In the present ex- 
amples, we just apply the third order polynomial ansatz 
with the boundary conditions in Eqs. (jT5)) - (|2"l"j) . so that 
the applied magnetic field Bq should not go beyond the 
upper limit in Fig. [2] As the upper limit field grows for 
smaller times tf, this is not a problem in practice. Figure 
[3] shows examples of spin flip for different values of B . 

With the functions 9 and <p fixed [see Fig. [3] (a)], the 
designed electric fields, £ x {t) and £ y (t), corresponding to 
B„ = 0.01 T and B = 0.07 T (close to the upper limit), 
are depicted in Fig. [3j(b) and (c). The populations (not 
shown) of the two spin states, given by Pi = cos 2 (9/2) 
and P-i = sin 2 (9/2), cross each other smoothly as 9 goes 
from to 7r. The choice of Bq determines the trajectory 
on the Bloch sphere for a given tf. When Bq approaches 
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FIG. 4: (Color online) Time evolution of cos# (a) and sirup 
(b) for the same Hamiltonian with the designed electric fields 
and Bo = 0.01 T, where e = 0, tp = tt/2 (solid blue); e = 0.01, 
<po = 7r/2 (dashed red); e = 0.01, v?o = 7r/4 (dot-dashed 
black). Other parameters are the same as in Fig. [3] 



the upper limit, the electric fields exhibit sharp peaks, see 
Fig. [3J (c). The smooth time-dependence in Fig. [3] (b) 
is well suited for the applications, while the complicated 
dependence in Fig. |3J (c) should be avoided. Undesirable 
excitation of the orbital modes does not occur here since 
the spin flip tf ~ 1 ns, while the energy split of the 
orbital states in typical QDs exceeds 0.1 meV. Therefore, 
regarding the orbital motion, our perturbation is strongly 
adiabatic, and no orbital excitation occurs. 

To check the stability with respect to initialization er- 
rors, we assume now the initial state as (VI — ee lipo , \/e) T 
with an arbitrary phase ipo and find 9 and ip from Eqs. 
(fTTjl and ([T2|) for the same designed electric fields (Fig. 
|4|. The final value 6(tf) depends on the error e, but in- 
sensitive to the initial phase (po, while <p{tf) is sensitive 
to both initial conditions. Since our goal is to realize the 
spin flip, the final if is irrelevant, and the experimental 
effort should focus on achieving a small error e. 

Decoherence and noise effects.- To show feasibility of 
our approach, we study the effects of noise and decoher- 
ence on the spin-flip fidelity We begin with a generic ap- 
proach for coupling to the incoherent environment, based 
on the conventional Lindblad formalism as can arise, e.g., 
from interaction with the conduction electron bath. The 
master equation reads (3ll | 



l .[H eS ,p]-l([a x , 



V"--.. 

tf = 1 ns \ 



t, = 0.1ns 



\ 



5 10 15 20 25 30 

FIG. 5: (Color online) Fidelity as a function of 7 for different 
times tf = 0.1 ns (solid red) and tf = 1 ns (dotted black). The 
fidelity estimated from perturbation theory is also compared, 
for tf =0.1 ns (dashed blue, undistinguished) and tf — 1 ns 
(dot-dashed orange). Bo = 0.01 T and other parameters are 
the same as in Fig. [2] 



where 7 is the dephasing rate. We introduce the Bloch 
vector with components u = pi-\ + p—11, v = — — 
P-u), and w = pn — and obtain from Eq. (|37[) 
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+ Wy, [Vy,P]\ + Wz, Wz,P]]) , 



(22) 



(23) 



We solve Eq. ([23| numerically and calculate fidelity 
F = \( — l|*(t/)>|, see Fig. \5\ For jt f < 1 the time- 
dependent perturbation theory [2l| yields the bound 
F > 1 — 27*/. Since the induced flip occurs very fast, it 
can overcome the main danger for the low-temperature 
spin manipulation in QDs coming from the hyperfine cou- 
pling to the nuclear spins, where the decoherence times 
exceed 100 ns 0- 

Another source of decoherence is the device-dependent 
noise in the electric field acting on the spin. This can 
be important when the relatively weak electric fields are 
applied. We analyze in detail the effect of this noise and 
find that our method is robust to this randomness in [32| . 

Conclusion and outlook.- We have proposed a fast and 
robust method to flip electron spin in a QD with SO 
coupling and weak perpendicular magnetic field. The 
spin-flip process, designed by Lewis-Ricsenfeld invari- 
ants, is faster than the decoherence for any known low- 
temperature dephasing mechanism. This method can be 
further complemented by optimal control theory for time- 
and energy-minimization subjected to different physical 
constraints j33| . Implementation of this technique will 
allow for high-fidelity spin-manipulation for quantum in- 
formation processing. 

We acknowledge funding by the Basque Government 
(Grants No. IT472-10 and BFI-2010-255), Ministerio 
de Ciencia e Innovacion (Grant No. FIS2009-12773- 
C02-01), UPV/EHU program (UFI 11/55), and Na- 
tional Natural Science Foundation of China (Grant No. 
61176118) and Shanghai Rising-Star Program (Grant No. 
12QH1400800). 
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I. Hamiltonian reduction we consider a 4-level system. The total Hamiltonian is 

H = Ho + H so + Hi nt . In the basis of two lowest spin-split 
Here we give the details on the derivation of the ef- orbital states, 
fective 2x2 Hamiltonian in the paper. For simplicity 



H 



/£i+A 2 /2 

E 1 - A z /2 
E 2 + A z /2 

\ E 2 - A z /2 j 



where A z is the Zeeman splitting, and 
/ 



^so + H nlt — 



-PA X 
'-a(iA x - A y ) 



'-a(iA x + A y ) (f3 ~ A x )p x - A y p y a(ip x + p y ) ^ 



: PA X 



<*(-iPx+Py) ~W + A x )p x - AyPy 



-pA a 
c 



(~P + A x )p x + A y p y -a(ip x + p y ) 
I a{ip x ~p y ) {P + A x )p x + A y py -a(iA x -A y ) 

I 



—a(iA x +A y ) 
c p 
-PA X 
c 



(24) 



(25) 



with p { = (ij>i\pi\ip2) and A, = eAt/mc. 

Our aim is to flip the spin between two levels taking 
into account the effects from orbital motion and the in- 
fluence of the other levels. As the energy gap of the 
orbital states is much larger than the Zeeman splitting 
A z , the Lowdin partition 28, [29j will enable us to ob- 
tain an effective 2x2 matrix Hamiltonian. We split the 
total Hamiltonian H into four blocks Q, B, C, C\ each 
of which is a 2 x 2 matrix, 



constant, the effective Hamiltonian becomes 



H 



Q C 
C*t B 



(26) 



The time-independent Schrddinger equation can be for- 
mally written as 



QF + CG = EF, 
C^F + BG = EG, 



(27) 
(28) 



Substituting the formal solution of Eq. (f2"5)l , G = (E — 
B)~ l C^F, into Eq. ([27]) . we get a closed equation for F, 



QF + C(E - BY^C^F = EF. 



(29) 



Therefore, the effective Hamiltonian for this subspace is 
given by 



H 



eff 



+ C(E - B)- 1 ^ 



including the effects from B, C and C'. 

Assuming that the electric field is not strong 
enough to excite other orbital states, namely, 
e max (\A x p x \ ,\Ayp y \) /mc <c E 2 — Ex, and keep- 
ing the first order term of the spin-orbit coupling 



H 0B = 



Ex + hlx + hxx h\ 2 

+ ^21 E\ 



hl2 



where 



hlx 



K 2 
hxi 



-h° 

n 22 



h 2 l = 
-h 2 2 

Ki = 



A z - -pA x , 
c 

e 

--a(iA x + A y ), 
c 

ej3 p x (A x p x + AyPy) 

mc E 2 — E\ 
eg (A x p x + AyPy) (ip x + p y ) 
mc E2 — Ei 



(30) 

(31) 
(32) 

(33) 

.(34) 



By a simple shift we may ignore the common term Ei on 
the diagonal, and finally express the Hamiltonian as 



2 



Z X + iY 
X-iY -Z 



(35) 



where X = B 2 (l + £„), Y = (a/0)(l + £ X )B U Z = 
B + (1 + L)Bi with Bi = -2epA x /(cgfx B ), B 2 = 
—2ej3A y /(cg^B)- According to perturbation theory, the 
effects from other states can be characterized by £ x and 
& = 2 En>i KV'ibilV'n}! 2 /[m(E n -Ex)], summing up 
all contributions. 



II. Apparatus Noise Effect 

Besides the decoherence resulting from the interaction 
between the system and the environment introduced in 
the main text, we consider apparatus noise, i.e. the 
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FIG. 6: (Color online) Dependence of fidelity on Ag, given by 
tf = 1 ns (solid black) and tf = 0.1 ns (dashed red). Other 
parameters are the same as in Fig. 2 (b) in the main text: 
ha = 2 x 10" 6 meV-cm, hfi = 10" 6 meV-cm, g = -0.44 for 
GaAs, Bo = 0.01 T, = 0. 



Hamiltonian H cS [Eq. ([35]) ] perturbed by a stochastic 



part H nolse describing noise from the electric field source. 
Therefore, the stochastic Schrodinger equation is 

ihj t ^(t) = (H cS + H noisc )y(t), (36) 

where H aoise = XH'C(t), (((t)) = 0, (C(*)C(*')) = *(*-f ) 
and A = Xo^/tJ, where Ao is the strength of noise. W 
is the part of H eS contributing from the time-dependent 
electric field. Then the density matrix [34| is 

P=-\\H cS , p]-^[H',[H',p}} (37) 

Following the same procedure to derive the Bloch equa- 
tion in the main text after Eq. (22), we obtain the corre- 
sponding equation with source noise. In order to analyze 
the effects from randomness induced by the electric field, 
we calculate numerically the fidelity F = |( — l|\l/(f^)|, 
see Fig. [6] The high fidelity shows that our method is 
robust to the source noise. 



